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^^ , Abstract. A relativistic model for quasielastic (QE) lepton-nucleus scattering is presented. The effects of final-state inter- 

Cn ■ actions (FSI) between the ejected nucleon and the residual nucleus are described in the relativistic Green's function (RGF) 

}_( ' model where FSI are consistently described with exclusive scattering using a complex optical potential. The results of the 

^1^' model are compared with experimental results of electron and neutrino scattering. 
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FINAL-STATE INTERACTIONS IN LEPTON-NUCLEUS SCATTERING 

Q I In the QE region the nuclear response to an electroweak probe is dominated by one-nucleon knockout processes, where 

^ ■ the scattering occurs with only one nucleon while the remaining nucleons of the target behave as simple spectators. 

Ph I The reaction can adequately be described in the relativistic impulse approximation (lA) by the sum of incoherent 

processes involving only one nucleon scattering and the components of the hadron tensor are obtained from the sum, 
over all the single-particle (s.p.) shell-model states, of the squared absolute value of the transition matrix elements 

^ . of the single-nucleon current. A reliable description of FSI is an essential ingredient for the comparison with data. 

^O ' The relevance of FSI has been clearly stated for the exclusive {e,e'p) reaction, where the use of complex optical 

^C , potentials (OP) in the distorted-wave impulse approximation (DWIA) is required [1,2, 3, 4, 5, 6, 7, 8]. The imaginary 

part of the OP produces an absorption that reduces the cross section and accounts for the fact that, if other channels 
are open besides the elastic one, part of the incident flux is lost in the elastically scattered beam and goes to the other 
^r i (inelastic) channels which are open. In the inclusive scattering only the emitted lepton is detected, the final nuclear 

state is not determined and all elastic and inelastic channels contribute. This requires a different treatment of FSI 
where all final-state channels should be retained and the total flux, although redistributed among all possible channels, 
must be conserved. Different approaches have been used to describe FSI in relativistic calculations for the inclusive 
QE electron- and neutrino-nucleus scattering. In the relativistic plane-wave impulse approximation (RPWIA), FSI are 
simply neglected. In another approach, FSI are accounted for in relativistic DWIA (RDWIA) calculations by including 

}J] , only the real part of the relativistic optical potential (rROP). 

C^_- In the RGF techniques [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19], FSI are described in flie inclusive scattering by the 

same complex OP as in the exclusive scattering, but the imaginary part is used in the two cases in a different way and 
in the inclusive reaction the flux, although is redistributed in all the channels, is conserved. In the RGF model with 
suitable approximations, which are mainly related to the impulse approximation, the components of the hadron tensor 
are written in terms of the s.p. optical model Green's function. The explicit calculation of the s.p. Green's function can 
be avoided by its spectral representation, which is based on a biorthogonal expansion in terms of the eigenfunctions of 
the non-Hermitian optical potential and of its Hermitian conjugate 

[^^T-r\E)]\xi'\E))^0, [^-T-y{E)]\x^-\E))=0 . (1) 

The expanded form for the s.p. expression of the hadron tensor components is [10, 1 1] 
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FIGURE 1. Differential cross section of the reactions '^0(e,e') for beam energy e — 1080 MeV and scattering angle j3 — 32° 
and Fe{e,e') for beam energy e = 2020 MeV and scattering angle d = 20°. Experimental data from [25] ( O) and [26] (^^Fe). 



where ^ denotes the principal value of the integral, n is the eigenstate of the residual nucleus with energy £„, and 
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and similar expressions for the terms with fJ. ^V. The factor •\/l — y'{E), where '^'{E) is the energy derivative of the 
optical potential, accounts for interference effects between different channels and allows the replacement of the mean 
field Y with the phenomenological OP. 

Disregarding the square root correction, the second matrix element in Eq. (3) is the transition amplitude for the 
single-nucleon knockout from a nucleus in the state | ^o) leaving the residual nucleus in the state | n) and it is similar 
to the usual DWIA expression for the transition amplitude of the exclusive single-nucleon knockout, i.e., the imaginary 
part of Y ' gives an attenuation of the strength that is inconsistent with the inclusive process, where all the inelastic 
channels must be considered and the total flux must be conserved. This compensation is performed by the first matrix 
element in the right hand side of Eq. (3), which involves the eigenfunction x^ '{E) of the Hermitian conjugate optical 
potential, where the imaginary part has an opposite sign and has the effect of increasing the strength. Therefore, in the 
RGF approach the imaginary part of the optical potential redistributes the flux lost in a channel in the other channels, 
and in the sum over n the total flux is conserved. The RGF model requires the calculations of matrix elements of the 
same type as in usual RDWIA models, but involves eigenfunctions of both Y and Y^: FSI are described by the 
same complex optical potential as in RDWIA thus providing a consistent treatment of FSI in the exclusive and in the 
inclusive scattering. 



RESULTS FOR ELECTRON AND NEUTRINO SCATTERING 



The first measurement of the charged-current quasielastic (CCQE) flux-averaged double-differential muon neutrino 
cross section on '^C in an energy range up to « 3 GeV that have been reported by the IVliniBooNE Collaboration 
[20] have raised extensive discussions. In particular, the experimental cross section is usually underestimated by the 
relativistic Fermi gas model and by other more sophisticated models based on the lA [21, 22, 23], unless the nucleon 
axial mass M^ is significantly enlarged with respect to the world average value of 1.03 GeV/c^. Despite the fact the 
larger axial mass obtained from the MiniBooNE data on carbon can be interpreted as an effective way to include 
medium effects which are not taken into account, it is clear that, before drawing conclusions, a precise knowledge of 
lepton-nucleus cross sections, where uncertainties on nuclear effects are reduced as much as possible, is necessary. 
Moreover, any reliable calculation for neutrino scattering should first be tested against electron scattering in the same 
kinematical conditions. 

In this section several results obtained with the RGF model for electron and neutrino-nucleus scattering are 
discussed. Some results of the RGF for inclusive electron and neutrino scattering are presented in [16, 18], where 
they are compared with the results obtained with the relativistic mean field model [24], which uses the same strong 
real potential already considered in describing the bound states to evaluate the scattering wave functions. 

As an example, in Fig. 1 the RGF results are compared with the experimental (e,e') cross section on ^^O and 
^^Fe. In all the calculations presented here the bound nucleon states are self-consistent Dirac-Hartree solutions derived 
within a relativistic mean-field approach [27], and in the RGF different parameterizations have been used for the 
relativistic optical potential: the energy-dependent and A-dependent (where A is the mass number) EDADl and the 
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FIGURE 2. Flux-averaged double-differential cross section per target nucleon for the CCQE '^C(v^,/i ) reaction calculated 
with the RGF-EDAI (dashed lines) and the RGF-EDADl (solid lines) displayed versus T^ for three bins of cos iJ^. The data are 
fromMiniBooNE[20]. 
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FIGURE 3. Total CCQE cross section per target nucleon as a function of Ey (left panel) and Ey (right panel) calculated with the 
RGF-EDADl (solid line), the RGF-EDAI (dashed hne), the rROP (dotted line), and the RPWIA (dot-dashed hne). The data are 
from MiniBooNE [20, 29]. 



energy-dependent but A-independent EDAI complex phenomenological potentials of [28]. The shape followed by 
the RGF cross sections fits well the slope shown by the data, in particular approaching the peak region, where the 
RGF produces cross sections in reasonable agreement with the data. Although satisfactory on general grounds, the 
comparison with data in Fig. 1 cannot be conclusive until contributions beyond the QE peak, like meson exchange 
currents and A effects, which may play a significant role in the analysis of data even at the maximum of the QE peak, 
are carefully evaluated. 

In Fig. 2 the CCQE double-differential ^^C{v^,jJ.^) cross sections averaged over the neutiino flux is displayed as 
a function of the muon kinetic energy T^ for three bins of cosi}^, where i>^ is the muon scattering angle. In all the 
calculations of neutrino-nucleus scattering the standard value of the nucleon axial mass, i.e., Mj^ = 1.03 GeV/c^ has 
been used. A good agreement with the MiniBooNE data of [20] is generally shown by the RGF cross sections [18]. 
The differences between the RGF results with the two optical potentials are enhanced in the peak region but they 
always are of the order of the experimental eiTors. The EDADl and EDAI potentials yield close predictions for the 
bin 0.4 < cos iJ^ < 0.5; small differences are seen in the bin 0.8 < cos d-^ < 0.9, while larger differences can be found 
for the bin 0.7 < cos d-^ < 0.8. Nevertheless, the RGF-EDAI cross section is always larger than the RGF-EDADl one. 
The differences between the RGF-EDAI and the RGF-EDADl results are due to the different imaginary parts. 

In Fig. 3 the total CCQE cross sections per target nucleon for neutrino and antineutrino scattering are displayed 
as functions of the neutrino or antineutrino energies Ey and Ey and compared with the "unfolded" MiniBooNE data 
[20, 29]. The RPWIA and rROP results usually underpredict the v^ data. Larger cross sections, in particular for larger 
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FIGURE 4. NCE flux-averaged cross section per target nucleon as a function of (^ calculated with the RGF-EDADl (solid line) 
and the RGF-EDAI (dashed hne). The data are from MiniBooNE [30]. 



values of Ey are obtained in the RGF with both optical potentials [18]. The differences between the RGF-EDAI and 
the RGF-EDAD 1 results, being RGF-EDAI in good agreement with the shape and magnitude of the v^ data, are due 
to the different imaginary parts. The enhancement of the RGF cross sections is due to the translation to the inclusive 
strength of the overall effect of inelastic channels which are not incorporated in other models based on the lA. 

The total CCQE cross section for V|j scattering is displayed in the right panel of Fig. 3. Also in this case the RGF 
results are usually larger than the RPWIA and rROP ones. The differences between the RGF-EDADl and RGF-EDAI 
results are significant but somewhat smaller than for neutrino scattering. A reasonable agreement with the experimental 
data is obtained, but the RGF-EDAI calculations are larger than the data up to £9 ~ 1 .5 GeV, while a better agreement 
is obtained with the RGF-EDADl ones. The different behavior of the cross sections for neutrino and antineutrino 
scattering is related to the relative strength of the vector-axial response, which is constructive in v scattering and 
destructive in v scattering, with respect to the longitudinal and transverse ones [19]. 

The MiniBooNE Collaboration has reported [30] also a measurement of the flux-averaged differential cross section 
as a function of the four-momentum transferred squared, Q^ = —q^q^L, for neutral-current elastic (NCE) neutrino 
scattering on CH2 in a Q^ range up to « 1 .65 (GeV/c)^. The analysis of v-nucleus NCE reactions introduces additional 
complications, as the final neutrino cannot be measured and a final hadron has to be detected: the cross sections are 
therefore semi-inclusive in the hadronic sector and inclusive in the leptonic one [31, 32, 33]. Different relativistic 
descriptions of FSI are compared with the NCE MiniBooNE data in [34], while in Fig. 4 we show our RGF results 
calculated with both EDAI and EDADl potentials. Also in this case, the RGF produces large cross sections that are 
in nice agreement with the data. However, we stress the RGF is appropriate for the inclusive scattering where only 
the final lepton is detected, and thus can take into account also contributions that are not included in the experimental 
cross sections. 



SUMMARY AND CONCLUSIONS 



A deep understanding of the neutrino-nucleus cross sections is very important for the determination of neutrino 
oscillation parameters. Reliable theoretical models are required where all nuclear effects are well under control. 
Within the QE kinematics domain, the treatment of FSI is an essential ingredient for the comparison with data. In 
this contribution the RGF model for the inclusive QE electron and neutrino-nucleus scattering has been discussed. 
This model was originally developed for QE electron scattering, successfully tested in comparison with electron- 
scattering data, and then extended to neutrino-nucleus scattering. In the RGF model FSI are described in the inclusive 
scattering by the same complex optical potential as in the exclusive scattering, but the imaginary part is used in the 
two cases in a different way and in the inclusive process it is responsible for the redistribution of the flux in all the 
channels and the conservation of the flux. The RGF model gives results that are usually larger than results of other 
models based on the impulse approximation and that are in fair agreement with the CCQE MiniBooNE cross sections 
without the need to increase the standard value of the nucleon axial mass. However, the use of phenomenological 
optical potentials, does not allow us to disentangle the role of different reaction processes and explain in detail the 



origin of the enhancement with respect of other models. The important role of contributions other than direct one- 
nucleon emission has been confirmed by different models; it has been observed that the neutrino-nucleus reaction 
at MiniBooNE can have significant contributions from effects beyond the lA in some kinematic regions where the 
experimental neutrino flux has significant strength. In Refs. [21, 35, 36, 37, 38, 39] the contribution of multinucleon 
excitations to CCQE scattering has been found sizable and able to bring the theory in agreement with the MiniBooNE 
cross sections without increasing the value of Ma- A critical review of nuclear effects in NCE and CCQE scattering 
is presented in [40]. Moreover, processes involving two-body currents, whose role is discussed in [41], should also be 
taken into account explicitly and consistently in a model to clarify the role of multinucleon emission. Fully relativistic 
microscopic calculations of two-particle-two-hole (2p-2h) contributions are extremely difficult and may be bound 
to model-dependent assumptions. The RGF results are also affected by uncertainties in the determination of the 
phenomenological optical potential. At present, lacking a phenomenological OP which exactly fulfills the dispersion 
relations in the whole energy region of interest, the RGF prediction is not univocally determined from the elastic 
phenomenology. A better determination of a phenomenological relativistic optical potential, which closely fulfills the 
dispersion relations, would be anyhow desirable and deserves further investigation. 
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